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Abstract 

Charged particle production in deep-inelastic ep scattering is measured with the HI 
detector at HERA. The kinematic range of the analysis covers low photon virtualities, 
5 < < 100 GeV^, and small values of Bjorken-x, 10"'^ < x < 10~^. The analysis 
is performed in the hadronic centre-of-mass system. The charged particle densities are 
measured as a function of pseudorapidity and transverse momentum (p^) in the range 
< ry* < 5 and < < 10 GeV differentially in x and Q^. The data are compared to 
predictions from different Monte Carlo generators implementing various options for hadro- 
nisation and parton evolutions. 



Submitted to Eur. Phys. J. C 



C. Alexa^, V. Andreev^^, A. Baghdasaryan^^, S. Baghdasaryan^^, W. Bartel^^, K. Begzsuren^^, 
A. Belousov^^, P. Belov^^, V. Boudry^^, 1. Bozovic-Jelisavcic^, G. Brandt^'', M. Brinkmann^^, 
V. Brisson^^, D. Britzger^\ A. Buniatyan^"', A. Bylinkin^'^''^^, L. Bystritskaya^"^, 

A.J. Campbell", K.B. Cantun Avila^^ F. Ceccopieri\ K. Cemy^i, V. Chekelian^e, 

J.G. Contreras^^, J. Cvach=^°, J.B. Dainton^^, K. Daum^^'"^^, E.A. De Wolf', C. Diaconu^i, 

M. Dobre^ V. Dodonov^^, A. Dossanov^^'^^, G. Eckerlin", S. Egli^^ E. Elsen", L. Favart^, 

A. Fedotov^^, R. Felst^^, J. Feltesse^", J. Ferencei^^, D.-J. Fischer^^, M. Fleischer^^, 

A. Fomenko^^, E. Gabathuler^*^, J. Gayler", S. Ghazaryan", A. Glazov", L. Goerlich^, 

N. Gogitidze2^ M. Gouzevitch"'^^ C. Grab^^, A. Grebenyuk", T. Greenshawl^ 

G. Grindhammer2^ S. Habib", D. Haidt", R.C.W. Henderson^^, E. Hennekemper^^ 
M. Herbst^^ G. He^■era2^ M. Hildebrandt=^^ K.H. Hiller^^ J. Hladky=^°, D. Hoffmann^i, 

R. Horisberger"'^, T. Hreus^, F Huber^^, M. Jacquet^^, X. Janssen'', L. Jonsson^°, H. Jung^^'^, 

M. Kapichine^ C. Kiesling^e, M. Kleinl^ C. Kleinwort", R. Kogler^^ R Kostka'^^ 

M. Kramer", J. Kretzschmar^*, K. Kriiger", M.RJ. Landon^^, W. Lange^*^, R Laycock^*^, 

A. Lebedev^s, S. Levonian", K. Lipka"■''^ B. List", J. List", B. Lobodzinski", 

R. Lopez-Fernandez^^, V. Lubimov^'^ ^ E. Malinovski^^, H.-U. Martyn\ S.J. Maxfield^^, 

A. Mehta^^, A.B. Meyer", H. Meyer^^, J. Meyer", S. Mikocki''', 1. Milcewicz-Mika''', 

A. Morozov^, J.V. Morris*^, K. MuUer"^", Th. Naumann^^, RR. Newman^, C. Niebuhr", 

D. Nikitin^ G. Nowak^ K. Nowakl2''^^ J.E. Olsson", D. Ozerov", R Pahl", V. Palichik^ 
M. Pandurovic^ C. Pascaud^^, G.D. Patel^^, E. Perez^°■^^ A. Petrukhin", I. Picuric^^, 

H. Pirumov^^ D. Pitzl", R. Placakyte"■^^ B. Pokomy=^\ R. Polifka^^'^^ V. Radescu"'^^ 
N. Raicevic^^, T. Ravdandorj"''', P. Reimer^°, E. Rizvi^^, P. Robmann'^°, R. Roosen'', 

A. Rostovtsev^'^, M. Rotaru^ J.E. Ruiz Tabasco2^ S. Rusakov2^ D. Salek3\ D.PC. Sankey^ 
M. Sauter^^ E. Sauvan^i-^s, S. Schmitt", L. Schoeffel^o, A. Schoningl^ 

H.-C. Schultz-Coulon^^ F. Sefkow", S. Shushkevich", Y. Soloviev"■2^ P Sopicki^ 
D. South", V. Spaskov'', A. Specka^^ Z. Staykova'', M. Steder", B. Stella^^ G. Stoicea^ 
U. Straumann''", T. Sykora^'^\ PD. Thompson^, D. Traynor^^ P Truol'^o, I. Tsakov^^ 

B. Tseepeldorj'^^'", J. Tumau^ A. Valkarova^i, C. Vallee^i, P Van Mechelen\ Y. Vazdik^^, 

D. Wegener^ E. Wunsch", J. Zacek^i, J. Zalesak^o, Z. Zhang^^, R. Zlebcfk^i, H. Zohrabyan^^, 
and F. Zomer^^ 



^ I. Physikalisches Institut der RWTH, Aachen, Germany 

^ Vinca Institute of Nuclear Sciences, University of Belgrade, 1 100 Belgrade, Serbia 

^ School of Physics and Astronomy, University of Birmingham, Birmingham, UK^ 

^ Inter- University Institute for High Energies ULB-VUB, Brussels and Universiteit Antwerpen, 

Antwerpen, Belgium'^ 

^ National Institute for Physics and Nuclear Engineering (NIPNE) , Bucharest, Romania'' 
^ STFC, Rutherford Appleton Laboratory, Didcot, Oxfordshire, UK^ 
^ Institute for Nuclear Physics, Cracow, Poland'^ 
^ Institut fUr Physik, TU Dortmund, Dortmund, Germany'^ 
^ Joint Institute for Nuclear Research, Dubna, Russia 
^° CEA, DSM/Irfu, CE-Saclay, Gif-sur-Yvette, France 
DESY, Hamburg, Germany 

Institut fur Experimentalphysik, Universitat Hamburg, Hamburg, Germany"^ 
Max-Planck-Institut fUr Kernphysik, Heidelberg, Germany 
Physikalisches Institut, Universitat Heidelberg, Heidelberg, Germany"' 



1 



Kirchhojf-Institut fur Physik, Universitdt Heidelberg, Heidelberg, Germany^ 

Institute of Experimental Physics, Slovak Academy of Sciences, Kosice, Slovak Republic^ 

Department of Physics, University of Lancaster, Lancaster, UK'' 

Department of Physics, University of Liverpool, Liverpool, UK'' 

School of Physics and Astronomy, Queen Mary, University of London, London, UK'' 

Physics Department, University of Lund, Lund, Sweden^ 

21 CPPM, Aix-Marseille Univ, CNRS/IN2P3, 13288 Marseille, France 

22 Departamento de Fisica Aplicada, CINVESTAV, Merida, Yucatan, Mexico^ 
2^ Departamento de Fisica, CINVESTAV IPN, Mexico City, Mexico^ 

2^ Institute for Theoretical and Experimental Physics, Moscow, Russia^ 

2^ Lebedev Physical Institute, Moscow, Russia 

2® Max-Planck-Institut fur Physik, Munchen, Germany 

2^ LAL, Universite Paris-Sud, CNRS/IN2P3, Orsay, France 

2^ LLR, Ecole Polytechnique, CNRS/IN2P3, Palaiseau, France 

2^ Faculty of Science, University of Montenegro, Podgorica, Montenegro' 

Institute of Physics, Academy of Sciences of the Czech Republic, Praha, Czech Republic^ 

^1 Faculty of Mathematics and Physics, Charles University, Praha, Czech Republic^ 

^2 Dipartimento di Fisica Universita di Roma Tre and INFN Roma 3, Roma, Italy 
Institute for Nuclear Research and Nuclear Energy, Sofia, Bulgaria 
Institute of Physics and Technology of the Mongolian Academy of Sciences, Ulaanbaatar, 

Mongolia 

Paul Scherrer Institut, Villigen, Switzerland 
Fachbereich C, Universitdt Wuppertal, Wuppertal, Germany 
Yerevan Physics Institute, Yerevan, Armenia 
DESY, Zeuthen, Germany 

Institut fur Teilchenphysik, ETH, Zurich, Switzerland^ 
'^^ Physik-Institut der Universitdt Zurich, Zurich, Switzerland^ 

^1 Also at Rechenzentrum, Universitdt Wuppertal, Wuppertal, Germany 

^2 Also at IPNL, Universite Claude Bernard Lyon 1, CNRS/IN2P3, Villeurbanne, France 

Also at CERN, Geneva, Switzerland 

Also at Ulaanbaatar University, Ulaanbaatar, Mongolia 

Supported by the Initiative and Networking Fund of the Helmholtz Association (HGF) under 
the contract VH-NG-401 and SO-072 

Also at Moscow Institute of Physics and Technology, Moscow, Russia 
'^'^ Also at Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7 

Also at LAPP, Universite de Savoie, CNRS/IN2P3, Annecy-le-Vieux, France 

Department of Physics, Oxford University, Oxford, UK'' 

^ Deceased 



" Supported by the BundesministeriumfUr Bildung und Forschung, FRG, under contract 
numbers 05H09GUF 05H09VHC, 05H09VHF 05H16PEA 

'' Supported by the UK Science and Technology Facilities Council, and formerly by the UK 
Particle Physics and Astronomy Research Council 

Supported by FNRS-FWO-Vlaanderen, IISN-IIKW and IWT and by Interuniversity Attraction 
Poles Programme, Belgian Science Policy 



2 



Partially Supported by Polish Ministry of Science and Higher Education, grant 
DPN/Nl 68/DESY/2009 
" Supported by VEGA SR grant no. 2/7062/27 
^ Supported by the Swedish Natural Science Research Council 

^ Supported by the Ministry of Education of the Czech Republic under the projects LC527, 
INGO-LA09042 and MSM002 1620859 

^ Supported by the Swiss National Science Foundation 
' Supported by CONACYT, Mexico, grant 48778-F 

^ Russian Foundation for Basic Research (RFBR), grant no 1329.2008.2 and Ros atom 

Supported by the Romanian National Authority for Scientific Research under the contract PN 
09370101 

' Partially Supported by Ministry of Science of Montenegro, no. 05-1/3-3352 



3 



1 Introduction 



Deep-inelastic scattering (DIS) processes at the ep collider HERA can access small values of 
Bjorken-x at low four momentum transfers squared Q"^ of a few GeV^. In the region of low x, 
characterised by high densities of gluons and sea quarks in the proton, the parton interaction 
with the virtual photon may originate from a cascade of partons emitted prior to the interaction 
as illustrated in figured] In perturbative Quantum Chromodynamics (QCD) such multi-parton 



Figure 1: Generic diagram for deep-inelastic ep scattering at small x. The transverse momenta 
of the emitted gluons are labeled as p^^i, while the proton longitudinal momentum fractions and 
the transverse momenta carried by the propagating gluons are denoted by and kT,i, respec- 
tively. 

emissions are described only within certain approximations valid in restricted phase space re- 
gions. At sufficiently large and not too small x the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) [IJ evolution equation is expected to be a good approximation. The DGLAP 
equation corresponds to a strong ordering of the transverse momenta of the propagator partons, 
, with respect to the proton direction, which implies strong ordering of the transverse mo- 
menta of the emitted partons, p^,* ^ PT,i+i, in the parton cascade from the proton towards 
the virtual photon. At small x the DGLAP approximation is expected to become inadequate 
and the Balitsky-Fadin-Kuraev-Lipatov (BFKL) [2] scheme may be more appropriate, which 
has no ordering in kx of the partons along the ladder. The Ciafaloni-Catani-Fiorani-Marchesini 
(CCFM) evolution (311 aims to unify the DGLAP and BFKL approaches. It introduces angular 
ordering of gluon emissions to implement coherence effects. At small x the CCFM evolution 
equation is almost equivalent to the BFKL approach, while it reproduces the DGLAP equations 
for sufficiently large x and Q"^. 

Measurements of the proton structure function F2{x, Q"^) pi] are well described by the Next- 
to-Leading-Order (NLO) or Next-to-Next-to-Leading Order (NNLO) DGLAP evolution [jSHHl, 
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suggesting that this observable may be too inclusive to exhibit signals for BFKL dynamics. De- 
viations from the kx ordering at HERA are observed in jet production (W. 10], transverse energy 
flow P111I12L forward jet production [iTBliTSl and measurements of forward tt^ production [[T6l . 
Studies of the transverse momentum spectrum of charged particles have been proposed in fH] 
as a more direct probe of the underlying parton dynamics. In that paper it has been shown with 
the help of QCD models that the high-p^ tail is sensitive to parton radiation, while the contri- 
bution from hadronisation is small. Previous measurements of the transverse momentum and 
pseudorapidity dependence of particle densities performed by the HI collaboration [181 were 
limited in statistical precision. 

This paper presents a study of charged particle production in ep collisions for 5 < < 
100 GeV^. The analysis is performed in the hadronic centre-of-mass frame, i.e. in the virtual 
photon-proton rest frame. The charged particle densities as a function of pseudorapidity and 
transverse momentum are measured differentially in x and Q^. 

The data set used for the analysis was collected with the HI detector in 2006 when positrons 
and protons collided with energies of 27.6 GeV and 920 GeV, respectively, corresponding to a 
centre-of-mass energy of ^/s = 319 GeV. The integrated luminosity of the data set is 88.6 pb^^, 
which is about seventy times larger than what was available for the previous HI measurement 
[flSil . This allows for a more detailed study of the dynamical features of parton evolution in the 
proton at small x. 



2 QCD models 

Parton cascade and hadronisation processes leading to charged particle production in ep colli- 
sions are modeled using different Monte Carlo (MC) programs . Brief descriptions of the MC 
event generators considered in this analysis are given below. 

• The Rapgap (T9\ MC generator matches first order QCD matrix elements to DGLAP 
based leading logarithm approximations for parton showers with strongly ordered trans- 
verse momenta of subsequently emitted partons. The factorisation and renormalisation 
scales are set to /x/ = /i^ = a/Q^ + Pr"^, where pT is the transverse momentum of the 
outgoing hard parton from the matrix element in the centre-of-mass frame of the hard 
subsystem. 

• The D JANGOH [|20l MC generator uses the Colour Dipole Model (CDM) as implemented 
in ARIADNE [[2T]|. which models first order QCD processes and creates dipoles between 
coloured partons. Gluon emission is treated as radiation from these dipoles, and new 
dipoles are formed from the emitted gluons from which further radiation is possible. 
The radiation pattern of the dipoles includes interference effects, thus modelling gluon 
coherence. The transverse momenta of the emitted partons are not ordered, producing a 
configuration similar to the BFKL treatment of parton evolution |[22l . 

• The Cascade [23J MC generator uses off-shell leading order QCD matrix elements, sup- 
plemented with parton emissions based on the CCFM evolution equation. The equation 
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requires an unintegrated gluon density (see [^3l), which takes the transverse momenta of 
the propagators into account. In contrast to the DGLAP evolution equation, the CCFM 
equation only contains gluon splitting g gg. 

• The Herwig++ [|241 MC program with the POWHEG (POsitive Weight Hardest Emission) 
option (25), combines the full matrix element including virtual corrections at 0{as) with 
a DGLAP-like parton shower simulation. The Herwig++ MC program uses the Coher- 
ent Parton Branching algorithm [|26l l27l which is based on colour coherence to suppress 
branchings outside an angular-ordered region of phase space. Here, final state radiation 
is angular ordered and initial state radiation is ordered in E ■ 9 ^ px, where E and 9 are 
the energy and polar angle of the radiated parton, respectively. 

• Photoproduction background is generated with the Phojet [|28l program, which uses a 
two-component dual parton model [|29l including diffractive processes and vector meson 
production. 

The Rapgap and Djangoh predictions are calculated using the CTEQ6L(L0) [|30l set of 
Parton Distribution Functions (PDF), while Cascade uses the default AO unintegrated gluon 
density set [[3T| . The predictions of Herwig-i-i- were obtained with the default PDF MRST 

02 NLO (3T\. To simulate hadronisation the Lund string fragmentation model [|33l is used, as 
implemented in JETSET [34J for Djangoh and Pythia [35J for both Rapgap and Cascade. 
The parameters of the Lund string fragmentation model used here are tuned to describe e+e^ 
results [|36l . The tuning was performed by the ALEPH collaboration using hadronic Z decay 
data and the Pythia6.1 simulation with Bose-Einstein correlations turned on. In addition, the 
tune obtained by the Professor tool [|37l using LEP data is also tested. Herwig-i-i- incorporates 
the cluster model [38] of hadronisation, in which colour-singlet clusters of partons form after 
the perturbative phase and then decay into the observed hadrons. 

Djangoh and Rapgap are also used together with the HI detector simulation in order 
to determine the acceptance and efficiency and to estimate the systematic uncertainties asso- 
ciated with the measurement. The programs are interfaced to HERACLES [|39l to simulate the 
QED-radiative effects. The generated events are passed through a detailed simulation of the 
HI detector response based on the Geant simulation program [40] and are processed using 
the same reconstruction and analysis program chain as for data. For the determination of the 
detector effects both the Rapgap and Djangoh predictions are studied. Both models describe 
all relevant control distributions reasonably well [41 J. To improve the determination of the de- 
tector corrections the transverse momentum and pseudorapidity of charged particles as well as 
inelasticity y, defined as y = Q"^ /{s ■ x), are reweighted to the data [41J. The reweighting is 
applied to the generator quantities. 

3 Experimental method 
3.1 HI detector 

A full description of the HI detector can be found elsewhere [|42] - l44l and only the components 
most relevant for this analysis are briefly mentioned here. The coordinate system of HI is 
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defined such that the positive z axis is pointing in the direction of the proton beam (forward 
direction) and the nominal interaction point is located at ^ = 0. The polar angle Q is then 
defined with respect to this axis. The pseudorapidity is defined to be = — ln(tan(^^/2)). 

Charged particles are measured within the central tracking detector (CTD) in the polar an- 
gle range 20° < 9 < 165°, which is also used to reconstruct the interaction vertex. The CTD 
comprises two large cylindrical jet chambers (CJCs), and the silicon vertex detector [I45ll46l . 
The CTD is operated inside a 1.16 T solenoidal magnetic field. The CJCs are separated by a 
cylindrical drift chamber which improves the z coordinate reconstruction. A cylindrical multi- 
wire proportional chamber WT\ . which is mainly used in the trigger, is situated inside the inner 
CJC. The trajectories of charged particles are measured with a transverse momentum resolution 
of <j{pt)/pt ~ 0.2%pr/GeV©1.5%. The forward tracking detector (FTD) [48] measures the 
tracks of charged particles at polar angles 6° < 6* < 25°. In the region of angular overlap, FTD 
and short CTD track segments are used to reconstruct combined tracks, extending the detector 
acceptance for well-reconstructed tracks. Both the CTD and the combined tracks are linked 
to hits in the vertex detectors: the central silicon tracker (CST) [i45ll46ll , the backward silicon 
tracker (BST) [49 1 and the forward silicon tracker (FST) [50 1. These detectors provide precise 
spatial coordinate measurements and therefore significantly improve the primary vertex spatial 
resolution. The CST consists of two layers of double-sided silicon strip detectors surrounding 
the beam pipe covering an angular range of 30° < 9 < 150° for tracks passing through both 
layers. The BST consists of six double wheels of strip detectors measuring the transverse coor- 
dinates of charged particles. The FST design is similar to the BST and consists of five double 
wheels of single- sided strip detectors. 

The lead- scintillating fibre calorimeter (SpaCal) (44] covering the region 153° < 9 < 177.5°, 
has electromagnetic and hadronic sections. The calorimeter is used to measure the scattered 
positron and the backward hadronic energy flow. The energy resolution for positrons in the 
electromagnetic section is a{E)/E ^ 7.1%/^y E/GeV © 1%, as determined in test beam mea- 
surements [51]. The SpaCal provides energy and time-of-flight information used for triggering 
purposes. A backward proportional chamber (BPC) in front of the SpaCal is used to improve 
the angular measurement of the scattered lepton. The liquid argon (LAr) calorimeter [52] covers 
the range 4° < < 154° and is used in this analysis in the reconstruction of the hadronic final 
state. It has an energy resolution of a{E)/E ^ 50%/ ^/E/GeV © 2% for hadronic showers, as 
obtained from test beam measurements [|53l . 

3.2 Event reconstruction 

The DIS kinematics is reconstructed based on the measurement of the scattered electron and the 
hadronic final state (HFS) particles. In the so-called eS-method [|54ll the kinematic variables 
Q^, y and x are given by: 



where s is the square of the centre-of-mass energy, E'^ and 9^ the energy and polar angle of 
the scattered lepton, respectively, E^ being the energy of incoming lepton and E = J^ii^i ~ 




and X 



(1) 



s ■ y 
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Pz^i) where the sum runs over all hadronic final state (HFS) particles i. This method provides 
an optimum in resolution of the kinematic variables and shows only little sensitivity to QED 
radiative effects. The HFS particles are reconstructed using an energy flow algorithm [55 1. This 
algorithm combines charged particle tracks and calorimetric energy clusters, taking into account 
their respective resolution and geometric overlap, into hadronic objects, while avoiding double 
counting of energy. 

3.3 Data selection 

DIS events were recorded using triggers based on electromagnetic energy deposits in the SpaCal 
calorimeter. The trigger efficiency is determined using independently triggered data. For DIS 
events the trigger inefficiency is negligible in the kinematic region of the analysis. 

The scattered lepton, defined by the most energetic SpaCal cluster, is required to have an 
energy E'^ larger than 12 GeV. The kinematical phase space is defined by 5 < < 100 GeV^ 
and 0.05 < y < 0.6, corresponding to the geometric acceptance of the SpaCal. The upper cut 
on y reduces background from photoproduction. In addition, x is required to be in the range of 
0.0001 <x < 0.01. 

Additional selections are made to reduce QED radiation effects and to suppress background 
events. The z coordinate of the event vertex is required to be within 35 cm of the nominal 
interaction point. Events with high energy initial state photon radiation are rejected by requiring 
35 < ~ Pz,i) < GeV. Here, the sum extends over all HFS particles and the scattered 

electron. This cut further suppresses photoproduction background events to a level of about 
0.5%. 

The tracks used in the analysis are measured in the CTD alone (central tracks) or result from 
combinations of CTD and FTD information (combined tracks). Central tracks are required to 
have transverse momenta in the laboratory frame pr > 150 MeV. The momentum of a combined 
track is required to be larger than 0.5 GeV to ensure that the track has enough momentum to 
cross the endwall of the CJC. Both central and combined tracks are required to originate from 
the primary event vertex and to be in the pseudorapidity range — 2 < < 2.5 measured in the 
laboratory frame. Using only tracks assigned to the event vertex, the contributions from in-flight 
decays of Kg, A and from photon conversions and from other secondary decays are reduced. 
Further track quality cuts [41] are applied to ensure a high purity of the track reconstruction. 

3.4 Definition of experimental observables 

The results of this analysis are presented in the hadronic centre-of-mass frame (HCM), to min- 
imise the effect of the transverse boost from the virtual photon. The transformation to the HCM 
frame is reconstructed with the knowledge of the kinematic variables and y ||4T1| . The trans- 
verse momentum and pseudorapidity of charged particles in the HCM frame are labelled as p^ 
and 7]*. Since in this frame the positive z* axis is defined by the direction of the virtual photon, 
HFS particles with r/* > belong to the current hemisphere and particles with r/* < originate 
from the target (proton remnant) hemisphere. 
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Charged particle densities as a function of transverse momentum and pseudorapidity are 
defined as (l/iV)(dn/dp^) and (l/iV)(dn/d?7*), respectively. Here, dn is the total number 
of charged particles with transverse momentum (pseudorapidity) in the dp^ (dr]*) bin and N 
denotes the number of selected DIS events. For distributions measured differentially in x and 
Q^, dn and N are the numbers for the respective (x, Q"^) bin. 

Hadronisation is expected to be more relevant at small transverse momenta, while the hard 
parton radiation is expected to contribute more significantly at high {p^ > 1 GeV) fTT}. To 
distinguish hadronisations effects from parton evolution signatures, the charged particle density 
is measured as a function of 77* for < < 1 GeV and for 1 < < 10 GeV. The p^ de- 
pendence of the charged particle densities is studied in two different pseudorapidity intervals, 
< 1]* < 1.5 and 1.5 < r]* < 5, referred to as the "central region" and "current region", re- 
spectively, as illustrated in figure |2] Such division approximately defines the regions where the 
sensitivity to the hard scatter is largest (current region) , and where the parton shower models 
can be tested (central region). The target region, 77* < 0, is not accessible in this analysis. 



P -direction 

<S 



Central Current 
< r|*< 1-5 1-5 < r|*< 5 

y- direction 




Figure 2: The two pseudorapidity regions analysed in this paper. The region < rj* < 1.5 and 
1.5 < ?7* < 5, are denoted as "central" and "current" regions, respectively. 



4 Data corrections 



The data are corrected to the number of stable charged particles including charged hyperons, 
with proper lifetime cr > 10 mm, in the phase space given in table [B Correction factors 
are calculated for each analysis bin from the ratio of the number of generated stable charged 
particles to the number of reconstructed tracks. The bin widths are chosen such that a purityQ 
of more than 75% is ensured in all bins. The correction takes into account detector effects 
like limited resolution and losses near the phase space boundaries, as well as a small residual 
contamination from weak decays of neutral particles (e.g. K'^ and A). 

In addition to migrations between bins inside the measurement phase space, there are migra- 
tions from outside of the analysis phase space and there is background from photoproduction. 

'The purity is defined as the ratio of the number of charged particles generated and reconstructed in a given bin 
to the total number of charged particles in the phase space of the analysis which are reconstructed in this bin. 
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DIS selection 


Four momentum transfers squared 


5 < g2 < 100 GeV^ 


Inelasticity 


0.05 <y <0.Q 


Charged particles 


Pseudorapidity in the laboratory frame 


-2<ri< 2.5 


Transverse momentum in the laboratory frame 


Pt > 150 MeV 


Pseudorapidity in the HCM frame 


<r]* <5 


Transverse momentum in the HCM frame 


< < 10 GeV 



Table 1: Phase space for charged particles. 



These contributions are subtracted prior to applying the correction factors according to the pro- 
cedure outlined in [[T5l . 

The Djangoh MC was used to correct the data. The differences to the correction factors 
obtained from Rapgap are taken as systematic uncertainties. The correction factors strongly 
depend on i]* and to a lesser extent on p^. In the 1.5 < r^* < 5 region they vary between 1 
and 1.8 with the largest values seen at high p^ and large i]*. In the < r^* < 1.5 region, the 
correction factors rise up to 2.6 at high p^, due to the limited detector acceptance in this region. 
The two MC models predict very similar correction factors for most of the phase space region, 
but differences up to 5.5% are observed at small r/* and high p^. 



5 Systematic uncertainties 

The following sources of systematic uncertainties are considered for all measured quantities. 

• The systematic uncertainty on the SpaCal energy scale is 1% fSG^, which results in a 
systematic uncertainty of typically 0.4% for the measured single differential distributions. 

• The SpaCal angular resolution of 1 mrad leads to a systematic uncertainty of about 0.1% 
for the measured distributions. 

• The hadronic energy scale uncertainty is known to a precision of 2% ll57l . The scale 
uncertainty enters into the uncertainty of the phase space calculations, which depend on 
E — Pz of the HFS, and also affects the boost to the HCM frame. The systematic effect 
on the present measurements is about 0.3%. 

• The systematic uncertainty arising from the model dependence of the data correction is 
taken as the difference of the correction factors calculated using Rapgap and Djangoh 
MC. The resulting uncertainty on the measurements varies between 0.2% and 5.5%. 

• The systematic uncertainty associated with the track reconstruction (e.g. track recon- 
struction efficiency, vertex reconstruction efficiency, weak decays and nuclear interaction 
uncertainties) is estimated to be: 
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- 1% for central tracks, determined from the analysis of curling tracks and from the 
analysis of secondary vertices of charged particles located in the material between 
the two CJCs and originating from interactions with the detector material. The 
nuclear interaction cross sections of pions and kaons is found to be smaller in the 
simulation than in data. After correcting for these deficits, the agreement in the track 
efficiency between data and MC is found to be better than 1%. 

- 10% for combined tracks [141 11481 . This was checked using all selected central tracks, 
as well as by using pions from K'g decays, as a function of transverse momentum 
and pseudorapidity. Consistent results are obtained from both samples showing 
agreement of data and MC within 10%. 

The systematic uncertainty associated with the track reconstruction is applied as an in- 
dependent uncorrelated uncertainty on every data point. The resulting uncertainty on the 
measurements varies between 1% and 5.4% and is 1.6% on average. An additional sys- 
tematic uncertainty of 0.2% is assigned due to the different K^g contamination seen in 
data and MC for both central and combined tracks. The corresponding effect arising from 
A contamination is expected to be negligible. 

• The systematic uncertainty on the remaining photoproduction background is estimated to 
be 30%. This results in an uncertainty on the measured densities up to 0.9% at small x and 
Q^. At large x and the contribution from photoproduction is small and its uncertainty 
is negligible. 

The systematic uncertainties shown in the figures and tables are calculated by adding all 
contributions in quadrature. The total systematic uncertainty for the single differential mea- 
surements is below 2.5% for most analysis bins. 



6 Results 

The measurements of the charged particle densities as a function of pseudorapidity and trans- 
verse momentum in the phase space summarised in table [Hare listed in tables |2] to |9] and shown 
in figures |3l to [TOl 

6.1 Charged particle densities as a function of pseudorapidity 

The charged particle densities as a function of rj* were measured separately for < 1 GeV 
and for 1 < < 10 GeV, as shown in figure |3l In the soft region, the pseudorapidity dis- 
tribution is almost flat in the 1-5 < r^* < 3 range with about 1.7 charged particles per unit of 
pseudorapidity. The distribution falls at small rj* due to the cut on pseudorapidity in the labo- 
ratory frame. In the hard region the distribution becomes more peaked near rj* = 2.5, with 
a maximum of 0.23 charged particles per unit of pseudorapidity. For 1 < < 10 GeV the 
density increases rather strongly up to 7]* ~ 2.5, a behaviour expected from the strong ordering 
of transverse momentum towards the hard scattering vertex. 



11 



Figure |3] also shows the predictions of the DGLAP-like model Rapg AP based on different 
PDF sets. In the soft region all NLO PDFs (HERAPDFl.O (51, CTEQ6.6 [59"], GRV98NLO 
[|60l ) show similar results and predict less particles compared to calculations using the default 
LO PDF set CTEQ6L(L0). All predictions are close to the data. At large p^, differences be- 
tween the NLO PDF sets are observed, with CTEQ6L(L0) being closest to the data, although 
the differences to the data are still larger than the differences between the various PDF predic- 
tions. Similar PDF uncertainties are observed when using the CDM model as implemented in 
Djangoh. 

To check the sensitivity to hadronisation effects, the Rapgap predictions obtained with 
three sets of fragmentation parameters are compared to the data in figured parameters tuned 
by ALEPH [36], by the Professor tuning tool [37J and defauk Pythia6.424 fragmentation 
parameters. Significant differences between these three samples are seen in the soft region, 
where the data are best described by the ALEPH tune. At large transverse momenta they give 
similar results but none of them describes the data. 

Predictions from models with different approaches for QCD radiation (see section 2) are 
shown in figure |5l The data are compared to the CDM model DjANGOH, the DGLAP-based 
MC Rapgap and Herwig-i-i- and the CCFM model Cascade. In the soft plj, region, Djan- 
goh and Rapgap describe the data within the PDF uncertainties (figure |3]). Herwig-i-i-, which 
uses the cluster fragmentation model, provides a reasonable description of the data in the cen- 
tral region. The effect of not using the POWHEG option in Herwig-i-i- also has been investi- 
gated. Only small differences were observed which are not considered further in this paper. 
Cascade predicts too high multiplicities in most of the measured rj* range. In the region of 
1 <p*rp <IQ GeV the best description of the data is achieved by Djangoh. Rapgap strongly 
undershoots the data in the central region. Herwig-i-+ predicts a spectrum which is even below 
the prediction of Rapgap. Cascade is significantly above the data in a wide range of rj* . 

The charged particle densities as a function of rj* are shown in figure |6]for <1 GeV in 
eight different intervals of x and Q"^ . The data are compared to predictions of the Djangoh, 
Rapgap, Herwig-i-i- and Cascade generators. Djangoh provides a good description of the 
data over the full kinematic range. In general the description of the data by RAPGAP is some- 
what worse, with overshooting the data by about 10% at low x. Herwig-i-i- predicts smaller 
charged particle densities than observed in data in most of the phase space with differences of 
the order of 10% at the highest Q^. CASCADE is significantly above the data for r/* < 3 in all 
(x, Q^) bins. 

In figure |7] the charged particle densities as a function of rj* are shown in (x, Q'^) intervals 
for 1 < < 10 GeV . The shape of the distributions changes with x and more strongly 
than what is observed for < 1 GeV (figure |6]). At small values of x and the measured 
distribution is less dependent on rj* compared to the region at high x and . None of the 
models describes all aspects of the data. In general Djangoh is closest to the data. However 
it fails to describe the data at low and medium x in the central pseudorapidity region, with 
downwards deviations of the order of 20%. The RAPGAP prediction is below the data, with the 
strongest deviation observed at small x and small rj* . Herwig-i-i- significantly undershoots the 
data. The prediction of CASCADE agrees reasonably well with the measurement at low x and 
Q"^, but overshoots the data significantly as x or Q'^ increases. 
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6.2 Charged particle densities as a function of transverse momentum 

In figured] the charged particle densities as a function of are shown for two pseudorapidity 
intervals, < r^* < 1.5 (central) and 1.5 < ?7* < 5 (current). The shapes of the measured 
distributions in the two pseudorapidity ranges are similar. The spectrum falls by more than 
four orders of magnitude from < 1 GeV to ~ 8 GeV. The measurements are compared 
to the predictions of the DJANGOH, Rapgap, Herwig++ and CASCADE generators. DJANGOH 
provides in general a good description of the data, while only at high in the current region 
deviations from the measurement are observed. The other models fail to describe the data, with 
the strongest deviations being observed in the central region. The ratio of RAPGAP to data 
shows a sharp drop at ~ 1 GeV. The spectra predicted by Herwig++ are even softer 
than those predicted by Rapgap. Cascade in general produces higher particle densities than 
measured. 

In figures [9] and [TOl the charged particle densities as a function of p^ are shown for eight 
(x, Q^) intervals for the central and the current region, respectively. In the central region the 
measurement shows a dependence on x, such that the number of soft particles is decreasing with 
increasing x for fixed Q"^. In the current region this effect is less pronounced. The Djangoh 
model provides in general a good description of the data over the full kinematic range in both 
pseudorapidity regions, degrading at high in the lowest (x, Q"^) bin. Significant deviations 
of the Rapgap predictions from data are observed in the central region at low x and Q"^. The 
description becomes somewhat better at larger values of x and Q"^. The same trend is observed 
for the current region, but the overall data description is better. Herwig++ fails to describe the 
measurements at high p^ in the whole phase space. At lowest x and the spectrum is much 
softer than the one obtained with RAPGAP, while at high x and both predictions are similar. 
Cascade describes the data in the lowest x and Q"^ bin at high only. 



7 Conclusion 

This paper presents a study of charged particle production in ep collisions at low Q'^ mea- 
sured with the HI detector. The kinematic range of the analysis covers low photon virtualities, 
h < < 100 GeV^, and small values of x, 10^"^ < x < 10^^. The analysis is performed in 
the hadronic centre-of-mass system. The charged particle densities as a function of pseudora- 
pidity {rj*) and transverse momentum {p^) are measured differentially in x and Q"^. The charged 
particle densities as a function of pseudorapidity show different shapes, depending on the p^ 
range. For < < 1 GeV the density of particles is approximately constant for 1 < rf* < 3.5, 
while for 1 < < 10 GeV the density increases with increasing rj* up to rj* ^ 2.5, a behaviour 
expected from the strong ordering of transverse momentum towards the hard scattering vertex. 
The charged particle densities as a function of transverse momentum show an x dependence at 
small T]* (0 < r]* < 1.5), such that the number of soft particles is decreasing with increasing x, 
while in the 1.5 < 77* < 5 range this effect is less visible. 

In order to relate the charged hadron spectra to the parton dynamics at small x, the data are 
compared to QCD models with different evolution approaches for simulating the parton cascade 
and with different hadronisation schemes. The data allow the validity of different models to be 
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tested. At small p^, the data are reasonably well described by Djangoh (based on the Colour 

Dipole Model), as well as by Rapgap (based on the DGLAP shower evolution). At high 
and at low r]*, Rapgap severely undershoots the data. The differences are most pronounced at 
lowest X and Q'^, and decrease with increasing x and Q'^ values. Herwig++ which is also based 
on DGLAP but uses a cluster fragmentation model is significantly below the data over the full 
phase space. Cascade (based on CCFM) gives a reasonable description only at the lowest x 
and Q^, but overall predicts higher charged particle densities than observed in data. The Colour 
Dipole Model implemented in D JANGOH is the best among the considered models and provides 
a reasonable description of the data. 
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Table 2: Charged particle densities as a function of 77* for < < 1 GeV with relative statis- 
tical (stat.) and systematic (sys.) uncertainties given in per cent. The phase space is defined in 
table 1. 
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Table 3: Charged particle densities as a function of rj* for 1 < < 10 GeV with relative sta- 
tistical (stat.) and systematic (sys.) uncertainties given in per cent. The phase space is defined 
in table 1. 
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Table 4: Charged particle densities as a function of t]* for < < 1 GeV for different and 
x intervals with relative statistical (stat.) and systematic (sys.) uncertainties given in per cent. 
The phase space is defined in table 1. 
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Table 4: continued 
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0.33 


3.4 


5 < g2 < 10 GeV2 


1.5-2.0 


0.1758 


0.31 


3.3 


0.00024 <x< 0.0005 


2.0-2.5 


0.1903 


0.30 


5.2 




2.5-3.0 


0.2037 


0.29 


1.9 




3.0-3.7 


0.1626 


0.38 


2.6 




3.7 - 5.0 


0.0311 


1.33 


8.0 




U.O — i.O 


n 1 o/i /I 


U.DD 


4.0 


5 < g2 < 10 GeV^ 


1.5-2.0 


0.1675 


0.32 


4.7 




on o 


n 1 fii 1 
U.ioii 


n Qi 
U.oi 


Q n 


0.0005 <x< 0.002 


2.5-3.0 


0.1686 


0.32 


2.3 




3.0-3.7 


0.0978 


0.51 


3.0 




3.7-5.0 


0.01071 


2.21 


4.7 




0.0 - 0.5 


0.1506 


0.42 


6.1 




0.5 - 1.0 


0.1764 


0.30 


2.1 


10 < < 20 GeV^ 




1.0 - 1.5 


0.1959 


0.28 


2.7 


0.0002 <x< 0.00052 


1.5-2.0 


0.2180 


0.26 


1.7 




2.0-2.5 


0.2444 


0.25 


1.9 




2.5-3.0 


0.2646 


0.23 


2.7 




3.0 - 3.7 


0.2334 


0.31 


1.4 




3.7-5.0 


0.0552 


1.03 


2.7 



Table 5: Charged particle densities as a function of ?]* for 1 < < 10 GeV for different Q'^ 
and X intervals with relative statistical (stat.) and systematic (sys.) uncertainties given in per 
cent. The phase space is defined in table 1. 
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Q , X intervals 


rf range 


IjN ■ dn/d77* 


Stat. 


sys. 










(%) 


(%) 




0.0- 


1.0 


0.123 


0.69 


8.3 


10 < < 20 GeV^ 


1.0- 

1.5 - 


1.5 
2.0 


0.1910 
0.2145 


0.30 
0.27 


3.4 
3.1 


0.00052 <x < 0.0011 


2.0 - 


2.5 


0.2409 


0.26 


3.6 




2.5 - 


3.0 


0.2415 


0.26 


2.6 




3.0- 


3.7 


0.1750 


0.37 


2.4 




3.7- 


5.0 


0.0315 


1.34 


3.7 




0.5- 


1.5 


0.1496 


0.58 


6.6 


10 < g2 < 20 GeV^ 


1.5 - 


2.0 


0.2086 


0.27 


3.8 




2.0- 


2.5 


0.2246 


0.26 


4.5 


0.0011 <x < 0.0037 


2.5 - 


3.0 


0.1997 


0.28 


3.8 




3.0 - 


3.7 


0.1051 


0.48 


3.8 




3.7 - 


5.0 


0.01028 


2.38 


7.7 




0.0- 


0.5 


0.1581 


0.42 


6.1 


20 < g2 < 100 GeV^ 


0.5- 
1.0- 


1.0 

1.5 


0.2105 

0.2479 


0.27 

0.23 


2.3 

2.5 


0.0004 <x< 0.0017 


1.5 - 
2.0 - 


2.0 
2.5 


0.2820 
0.3188 


0.22 
0.20 


2.1 
2.3 




2.5 - 


3.0 


0.3386 


0.20 


1.4 




3.0- 


3.7 


0.2601 


0.28 


1.5 




3.7- 


5.0 


0.0602 


0.95 


1.8 




0.0- 


1.0 


0.118 


0.63 


9.0 


20 < g2 < 100 GeV2 


1.0- 
1.5- 


1.5 
2.0 


0.251 
0.2966 


0.22 
0.19 


4.5 
3.4 


0.0017 < X < 0.01 


2.0- 


2.5 


0.3167 


0.19 


2.1 




2.5- 


3.0 


0.2799 


0.20 


2.1 




3.0- 


3.7 


0.1587 


0.34 


2.4 




3.7- 


5.0 


0.02112 


1.49 


4.6 



Table 5: continued 
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< 77* < 1.5 


range 


1/N ■ dn/dp*j. 


Stat. 


sys. 


[GeV] 


[GeV-^] 


(%) 


(%) 


0.2 - 0.4 


3.952 


0.01 


2.0 


0.4 - 0.6 


2.431 


0.02 


1.6 


0.6-1.0 


0.954 


0.04 


1.8 


1.0-2.0 


0.1686 


0.15 


2.5 


2.0-4.0 


1.549 • 10-2 


0.70 


2.0 


4.0 - 10.0 


7.15 • 10-^ 


5.41 


1.9 



Table 6: Charged particle densities as a function as a function of in the region < ?7* < 1.5 
shown with relative statistical (stat.) and systematic (sys.) uncertainties given in per cent. The 
phase space is defined in table 1. 



1.5 < 77* < 5 


p^ range 


1/N ■ dn/dpT 


Stat. 


sys. 


[GeV] 


[GeV-i] 


(%) 


(%) 


0.0-0.3 


5.24 


0.01 


1.8 


0.3-0.6 


6.10 


0.01 


1.7 


0.6-1.0 


2.234 


0.02 


1.8 


1.0-1.5 


0.6193 


0.05 


1.5 


1.5 - 2.1 


0.1849 


0.10 


1.5 


2.1-3.0 


5.23 • 10-2 


0.23 


2.0 


3.0-4.0 


1.381 • 10-2 


0.47 


2.0 


4.0-5.0 


4.14 • 10-=^ 


0.84 


2.4 


5.0-6.3 


1.402 ■ 10-3 


1.67 


2.8 


6.3-7.9 


3.98- IQ-^ 


3.47 


2.5 


7.9 - 10.0 


1.061 • 10-^ 


7.60 


3.2 



Table 7: Charged particle densities as a function as a function of p^ in the region 1.5 < 77* < 5 
shown with relative statistical (stat.) and systematic (sys.) uncertainties given in per cent. The 
phase space is defined in table 1. 
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< T7* < 1.5 


Q^, X intervals 


range 


1/N ■ dn/dp^ 


Stat. 


sys. 




[GeV] 


[GeV-i] 


(%) 


(%) 




0.2 - 0.4 


4.76 


0.03 


3.2 


< < iU LreV 


n /I n 
U.4 — U.O 


o no 


U.U4 


z.y 




0.6-1.0 


1.144 


0.10 


3.7 


0.0001 <x < 0.00024 


1.0 - 2.0 


0.1955 


0.39 


3.0 




2.0-4.0 


1.489 • 10-2 


1.89 


3.2 




4.0-10.0 


5.69 • 10-^ 


15.36 


6.0 




0.2 - 0.4 


3.99 


0.05 


2.5 


5 < g2 < 10 GeV^ 


0.4 - 0.6 


2.53 


0.06 


2.6 




0.6-1.0 


0.994 


0.14 


2.9 


0.00024 <x< 0.0005 


1.0-2.0 


0.1611 


0.52 


3.3 




2.0-4.0 


1.286 • 10-2 


2.58 


2.7 




4.0-10.0 


5.1 • 10-^ 


21.40 


4.9 




0.2-0.6 


2.097 


0.11 


3.0 


5 < g2 < 10 GeV^ 


0.6-1.0 


0.659 


0.19 


3.4 


0.0005 <x< 0.002 


1.0-2.0 


0.1113 


0.74 


2.8 








2.0-4.0 


8.82 ■ 10-^ 


3.40 


3.1 




4.0-10.0 


3.3 • 10-^ 


34.47 


5.5 



Table 8: Charged particle densities as a function of p^ in the region < 77* < 1.5 for different 
and x intervals shown with relative statistical (stat.) and systematic (sys.) uncertainties 
given in per cent. The phase space is defined in table 1. 
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Q^, X intervals 


range 


l/N ■ dn/dpT 


Stat. 


sys. 




[GeV] 


[GeV-i] 


(%) 


(%) 




0.2-0.4 


4.97 


0.04 


2.5 


10 < g2 < 20 GeV2 


0.4-0.6 


3.060 


0.05 


2.7 




0.6-1.0 


1.229 


0.11 


2.6 


0.0002 <x< 0.00052 


1.0-2.0 


0.2155 


0.40 


3.5 




2.0-4.0 


1.960 • 10-2 


1.87 


3.1 




4.0-10.0 


9.2 • 10-^ 


14.13 


3.3 




0.2 - 0.4 




0.05 


o.u 


10 < g2 < 20 GeV^ 


0.4 - 0.6 


2.593 


0.06 


2.9 




0.6-1.0 


1.033 


0.13 


3.5 


0.00052 <x< 0.0011 


1.0-2.0 


0.1811 


0.50 


4.2 




2.0-4.0 


1.623 ■ 10-2 


2.34 


5.8 




4.0-10.0 


7.5 ■ 10-^ 


18.91 


3.2 




0.2-0.6 


2.124 


0.10 


4.5 


10 < < 20 GeV2 


0.6 - 1.0 


0.692 


0.18 


4.1 


0.0011 < X < 0.0037 


1.0-2.0 


0.129 


0.66 


4.8 








2.0-4.0 


1.119 ■ 10-2 


3.12 


3.0 




4.0-10.0 


4.42 • 10-^ 


28.92 


6.0 




0.2 - 0.4 


5.00 


0.04 


2.5 


20 < g2 < 100 GeV^ 


0.4 - 0.6 


3.156 


0.04 


2.4 




0.6-1.0 


1.296 


0.97 


2.8 


0.0004 <x< 0.0017 


1.0-2.0 


0.2474 


0.36 


3.7 




2.0-4.0 


2.579 • 10-2 


1.57 


4.7 




4.0-10.0 


1.40 • 10-=^ 


11.28 


2.8 




0.2-0.4 


3.356 


0.05 


2.8 


20 < g2 < 100 GeV2 


0.4-0.6 


2.254 


0.05 


2.8 




0.6 - 1.0 


0.945 


0.12 


3.6 


0.0017 <x< 0.01 


1.0-2.0 


0.1905 


0.42 


5.4 




2.0-4.0 


2.185 ■ 10-2 


1.74 


6.9 




4.0-10.0 


1.21 ■ 10-3 


13.18 


3.5 



Table 8: continued 
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1.5 < r/* < 5 


, X intervals 


range 


yN ■ dn/dp^ 


Stat. 


sys. 




rGeVl 


rGeV"M 


(%) 


(%) 




0.0 - 0.3 


5.60 


0.04 


2.5 




0.3 - 0.6 


6.83 


0.04 


2.8 




0.6 - 1.0 


2.290 


0.07 


3.6 


5 < Q2 < 10 GeV^ 


1.0 - 1.5 


0.639 


0.14 


3.3 


0.0001 <x< 0.00024 


1.5 - 2.1 


0.1897 


0.30 


2.7 




2.1 - 3.0 


5.12 • 10-2 


0.69 


4.0 




3.0 - 4.0 


1.378 • 10-2 


1.40 


4.8 




4.0-5.0 


3.97 • 10-3 


2.46 


8.6 




5.0 - 6.3 


1.40 • 10-3 


4.68 


9.7 




6.3 - 7.9 


4.34 • 10-^ 


9.33 


6.2 




7.9 - 10.0 


9.3 • 10-^ 


20.80 


13.2 




0.0 - 0.3 


5.57 


0.04 


2.8 




0.3 - 0.6 


6.64 


0.04 


2.6 




0.6 - 1.0 


2.207 


0.08 


2.2 


5 < Q2 < 10 GeV^ 


1.0 - 1.5 


0.587 


0.17 


1.7 


0.00024 < a; < 0.0005 


1.5 - 2.1 


0.161 


0.36 


1.8 




2.1 - 3.0 


4.12 ■ 10-2 


0.85 


3.2 




3.0 - 4.0 


1.002 ■ 10-2 


1.83 


5.7 




4.0-5.0 


2.63 • 10-3 


3.46 


6.7 




5.0 - 6.3 


7.69 • 10-^ 


7.08 


5.7 




6.3 - 7.9 


2.46 • 10-^ 


13.98 


10.4 




7.9 - 10.0 


6.65 • 10-^ 


32.90 


18.6 




0.0 - 0.3 


4.96 


0.05 


2.2 




0.3 - 0.6 


5.74 


0.04 


2.1 




0.6 - 1.0 


1.863 


0.08 


3.4 


5 < g2 < 10 GeV2 


1.0-1.5 


0.472 


0.19 


2.4 


0.0005 <x< 0.002 


1.5 - 2.1 


0.1214 


0.40 


2.5 




2.1-3.0 


2.85 ■ 10-2 


1.02 


5.5 




3.0-4.0 


6.14 ■ 10-3 


2.40 


6.1 




4.0-5.0 


1.55 ■ 10-3 


4.94 


5.4 




5.0-6.3 


4.401 ■ 10-^ 


10.49 


12.8 




6.3-7.9 


1.11 ■ 10-^ 


24.17 


13.6 



Table 9: Charged particle densities as a function of in the region 1.5 < r^* < 5 for different 
Q2 and X intervals shown with relative statistical (stat.) and systematic (sys.) uncertainties 
given in per cent. The phase space is defined in table 1. 



27 



Q"^, X intervals 


range 


V^V ■ dn/dp^ 


Stat. 


sys. 




[GeV] 


[GeV-^] 


(%) 


(%) 




0.0-0.3 


5.975 


0.04 


1.5 




0.3-0.6 


7.16 


0.04 


2.0 




0.6-1.0 


2.517 


0.07 


1.7 


10 < Q2 < 20 GeV2 


1.0-1.5 


0.742 


0.15 


2.4 


0.0002 <x < 0.00052 


1.5 - 2.1 


0.2317 


0.29 


2.3 




2.1 - 3.0 


6.82 ■ 10-2 


0.64 


3.4 




3.0-4.0 


1.848 ■ 10-2 


1.30 


3.6 




/in n 


O.OO ■ iU 


9 97 
Z.Z 1 


K 9 




5.0 - 6.3 


2.00 ■ 10-3 


4.45 


2.5 




6.3-7.9 


5.88 • 10-^ 


8.91 


2.9 




7.9 - 10.0 


1.942 • 10-^ 


17.77 


5.3 




0.0-0.3 


5.824 


0.05 


1.5 




0.3-0.6 


6.77 


0.04 


1.7 




0.6-1.0 


2.332 


0.08 


2.2 


10 < g2 < 20 GeV2 


1.0-1.5 


0.657 


0.16 


1.8 


0.00052 <x < 0.0011 


1.5 - 2.1 


0.1948 


0.31 


2.3 




2.1 - 3.0 


0.0538 


0.74 


3.0 




3.0-4.0 


1.297-10-2 


1.56 


5.9 




4.0 - 5.0 


3.94 ■ 10-3 


2.90 


4.5 




5.0 - 6.3 


1.33 • 10-3 


5.80 


6.9 




6.3 - 7.9 


3.40 ■ 10-"^ 


12.91 


6.3 




7.9 - 10.0 


7.7-10-5 


30.96 


8.6 




0.0-0.3 


5.12 


0.05 


2.2 




0.3-0.6 


5.65 


0.04 


2.4 


10 < (5^ < 20 GeV^ 


0.6-1.0 


1.946 


0.08 


2.0 




1.0 - 1.5 


0.538 


0.17 


2.2 


0.0011 <x< 0.0037 


1.5 - 2.1 


0.1481 


0.35 


3.9 




2.1-3.0 


3.89 • 10-2 


0.85 


3.6 




3.0-4.0 


8.643 • 10-3 


1.94 


3.4 




4.0-5.0 


2.10 • 10-3 


3.93 


4.7 




5.0-6.3 


6.43 ■ 10-^ 


8.90 


12.9 




6.3-10.0 


7.7-10-5 


46.70 


10.6 



Table 9: continued 
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(5^, X intervals 


range 


1/N ■ dn/dp*j. 


Stat. 


sys. 




[GeV] 


[GeV"^] 


(%) 


(%) 




0.0 - 0.3 


6.25 


0.04 


1.6 




0.3 - 0.6 


7.32 


0.04 


1.7 




n in 
u.o — i .u 




U.U ( 


9 1 
Z.L 


20 < g2 < 100 GeV^ 




1.0 - 1.5 


0.851 


0.13 


2.5 


0.0004 <x < 0.0017 


1.5 - 2.1 


0.2898 


0.24 


2.4 




2.1 - 3.0 


9.21 ■ 10^2 


0.52 


2.9 




3.0-4.0 


2.790 • 10-2 


1.01 


2.2 




4.0-5.0 


9.30 • 10-=^ 


1.75 


3.0 




5.0-6.3 


3.28 • 10-3 


3.37 


3.4 




6.3-7.9 


9.64 • 10-^ 


6.79 


4.2 




7.9 - 10.0 


3.06 • 10 * 


14.41 


3.9 




0.0-0.3 


5.652 


0.04 


1.5 




0.3-0.6 


6.052 


0.04 


1.6 




0.6- 1.0 


2.244 


0.06 


1.7 


20 < < 100 GeV2 


1.0- 1.5 


0.697 


0.13 


1.8 


0.0017 < X < 0.01 


1.5 - 2.1 


0.2323 


0.24 


1.9 




2.1 - 3.0 


7.12 ■ 10-2 


0.52 


2.1 




3.0-4.0 


2.025 ■ 10^2 


1.05 


1.9 




4.0-5.0 


6.32 • 10-3 


1.92 


3.9 




5.0-6.3 


2.20 • 10-3 


3.87 


4.0 




6.3-7.9 


5.91 • 10-^ 


8.42 


6.0 




7.9 - 10.0 


1.38 • 10-^ 


19.37 


20.1 



Table 9: continued 
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Figure 3: Charged particle density as a function of rj* for (a) < 1 GeV and for (b) 
1 < < 10 GeV compared to RAPGAP predictions with different proton PDFs. The predic- 
tions are obtained using the ALEPH tune. 
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Figure 4: Charged particle density as a function of rj* for (a) < 1 GeV for (b) 
1 < < 10 GeV compared to RAPGAP predictions for three different sets of fragmentation 
parameters. The predictions are obtained using CTEQ6L(L0) PDF. 
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Figure 5: Charged particle density as a function of rj* for (a) < 1 GeV for (b) 
1 < p:^ < 10 GeV compared to Djangoh, Rapgap, Herwig++ and CASCADE Monte Carlo 
predictions. 
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Figure 6: Charged particle density as a function of r]* for < 1 GeV for eight intervals of 
and X compared to Djangoh, Rapgap, Herwig++ and Cascade Monte Carlo predictions. 
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Figure 7: Charged particle density as a function of rj* for 1 < < 10 GeV for eight intervals of 
Q'^ and x compared to D jangoh, Rapgap, Herwig++ and Cascade Monte Carlo predictions. 



33 



(a) 



H1 data 

DJANGOH 

RAPGAP 



10 

> 

a 1 



CASCADE 
Herwig++ 



c 
■a 



I- 

10"^ 
10-^ 



^ 1.5 

Q 

0.5 



< Tj* < 1 .5 



H1 



I .... I .... I .... I . 



(b) 



H1 data 

DJANGOH 

RAPGAP 



CASCADE 
Herwig++ 



0123456789 1( 

P/ [GeV] 




123456789 1( 

P/ [GeV] 



Figure 8: Charged particle density as a function of in the ranges (a) < r^* < 1.5 and 
(b) 1.5 < rj* < 5 compared to DJANGOH, RAPGAP, Herwig++ and Cascade Monte Carlo 
predictions. The ratios of MC predictions to the measurements are shown on the bottom of the 
figure. 



34 



10 

1 

10-1 

10-2 

10-^ 

lO"* 
10 

1 

10-1 
10-2 
10-^ 

lo-'* 

10 

1 

10-1 

10-2 

10-^ 
10-" 



5 <Q <10GeV' 
0.0001 < X < 0.00024 



—•—I 

^ 



_±_ 



_±_ 



_±_ 



5 <Q <10GeV' 
0.00024 <x < 0.0005 



_±_ 



_±_ 



_±_ 



5 <Q <10GeV' 
0.0005 < X < 0.002 



HI 



_i_ 



10<Q <20 GeV^ 
0.0002 < X < 0.00052 



0.00052 <x < 0.0011 



10<Q <20 GeV^ 
0.001 1 < X < 0.0037 



_i_ 



_i_ 



_i_ 



_i_ 



_i_ 



_i_ 



20 <Q' <100GeV' 
0.0004 <x < 0.0017 



20 <Q^ <100GeV' 
0.0017 <x <0.01 



2 



8 2 



8 



< ri* < 1 .5 

>- HI data 
--■ DJANGOH 
RAPGAP 

— Herwig++ 

- - CASCADE 

1 ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ 

2 4 6 8 10 

P/ [GeV] 



Figure 9: Charged particle density as a function of in the range < r^* < 1.5 for eight 
intervals of Q"^ and x compared to Djangoh, Rapgap, Herwig++ and Cascade Monte Carlo 
predictions. 
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Figure 10: Charged particle density as a function of in the range 1.5 < rj* < 5 for eight 
intervals of Q"^ and x compared to Djangoh, Rapgap, Herwig++ and Cascade Monte Carlo 
predictions. 
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